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ABSTRACT 

Production of heavy quark antiquark systems in high energy heavy ion 
collisions must involve relativistic momentum components in a quantum me- 
chanical approach. If the color forces are screened in a deconfining medium, 
one can define the analog of a formation or separation time by an overlap 
integral in the nonrelativistic bound state rest frame. This time parameter 
has some interesting properties which depend on the momentum spectrum of 
the initial quarks. Consequences of these properties for the phenomenology of 
deconfinement signals are discussed. 

1. Introduction 

The motivation for this study is related to the possibility of using suppression 
of heavy quarkonia states in high energy heavy ion collisions as a signature for the 
formation of a quark-gluon plasmas The heavy quark-antiquark pair is produced by 
hard collisions of partons during the initial interaction times in the reactions. If a quark- 
gluon plasma is formed, one would expect that the confining color forces would be 
screened during the plasma lifetime t p and hence the heavy quark-antiquark pair could 
separate to a relative distance greater than that of the ordinary bound states in the 
confining potential. In the literature the minimum required time has been associated 
wiht the formation time of the bound state, but in this context it is more properly 
understood as a separation time t s in a nonconfining screened potential. When the 
confining potential reappears, the quark-antiquark pair are separated too far to fit 
into the confining potential region, so that at hadronization they are most likely to 
recombine with ordinary light quarks, thus leading to a suppression of the bound states 
with hidden flavor content and an enhancement of the open flavor states. 

2. Classical Interpretation 

The observation of this predicted suppression in J/ip production in O-U inter- 
actions by the NA38 collaboration at CERNcI has lead to an avalanche of theoretical 
and phenomenological papers. One feature of the data was immediately recognized 
as significant. The suppression was maximal for small J/ip transverse momentum P t 
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and gradually disappeared at a critical value P c ~ Mj/^ . This is immediately under- 
stood in simple terms, since the quark- ant iquark system can avoid suppression if its 
separation time (Lorentz dilated in the lab frame) is greater than either the plasma 
lifetime t p or the time of transit of the pair to the spatial boundary of the plasma 
region tf, = xE/P t where E is the transverse energy and x the transverse distance to 
the boundary. One must average over the production position of the quark-antiquark 
pair, which leads directly to a linear increase in the suppression factor. It reaches unity 
at P c = with d the transverse size of the plasma. If the plasma lifetime is the 
limiting parameter, the linear rise will be truncated by an immediate saturation at a 



smaller P c = MJ% - 1. Of course, the discontinuous values and slopes in this picture 

n 

are artifacts of the simple one-dimensional model. Realistic calculations™ found good 
agreement with the data for spatial parameters determined by the collision geometry, 
if the plasma lifetime t p ph lfm and the separation time t s ~ 0.7 fm. 

3. Quantum Mechanical Scenario 

These results seem to place quite severe constraints on the parameters of a pos- 
sible quark-gluon plasma. However, the dynamics of formation of the quark-antiquark 
pair in the hard collision tell a different story. The dominant mechanism for the pro- 
duction is gluon-gluon fusion, but any such process is characterized by a scale set by 
the heavy quark or heavy bound state mass. Hence the spacetime region involved in 
the production is Ax ~ -^p ~ Since most of the events are produced in a region 

V < m J/i> ~ 2mQ, one has > 1 > i-e. we are in the quantum-mechanical regime where 
there are not well-defined classical trajectories, and one cannot rely on parameters ex- 
tactly specifiying space and time events. One notes in addition that since the mass scale 
is set by the bound state or quark masses, one must also use relativistic momentum 
components in the quantum-mechanical wave packet which describes the position of 
the quarks. The procedure is then straightforward: one replaces the classical trajecto- 
ries of the quark and antiquark by expanding wave packets which are initially localized 
in space and propagate with central momenta provided in the production process. The 
probability of this pair forming the bound J/ip is just the overlap of the product of the 
wave packets onto a superposition of bound state times total momentum eigenstates. 
If the wave packet propagation is in the normal confining potential, these probabilities 
are time independent and the "formation time" plays no role. However, if the color 
forces are screened away as in a quark-gluon plasma, the evolution of the wave packets 
is altered and the probability of production of the J/tp becomes time-dependent. This 
is where the quantum mechanical analog of a separation time enters. 

4. A Simple Model 

We consider a one-dimensional example which can be done analytically. The 
initial wave packets for quark and antiquark are Gaussian with width a and initial 
(central) momenta p^. The bound state is also Gaussian with width erg and momentum 
P. The packets are allowed to propagate with mass m in a region of zero potential 
(simulating the plasma phase) for a time t, and then the probability amplitude for the 



formation of the bound state ap(t) is calculated: 
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The suppression factor is the ratio of that probability to its value in the absence 
of the plasma phase. One can define the quantum-mechanical separation time t s as the 
time required for this ratio to decrease by a factor of 2, yielding (for a non-relativistic 
bound state) 

t s = v / 3m(2a 2 + a 2 B ) (2) 

Alternatively, one can use the time-dependent shape of this function directly in 
the classical formulas for boundary crossing or plasma lifetime constraints, suitably 
Lorentz dilated for P^O. 

This procedure has been followed both in this simple model and for a full three- 
dimensional case with realistic collision geometry p and results in an effective classical 
separation time t s ® M ~ 0.2 — OAfm. According to the relationship between the effective 
parameters and transverse-momentum cutoff values, this then predicts P c values at least 
a factor of 2 or 3 greater than allowed by experiment. In retrospect, this is perhaps not 
so surprising, since the larger values of classical t s were estimated by average separation 
momenta between quark and antiquark in the bound state potential, and here we are 
allowing them to separate freely in zero relative potential as would be the case for a 
completely screened confining force. 

One must extend this analysis one more step, to study the situation when the 
total momentum of the bound state is relativistic. Up to now we have calculated the 
overlap integrals in the bound state rest frame, but if it is moving relativistically in 
the lab frame where the wave packet widths are defined, then one would expect that 
some Lorentz contraction of the widths should be included. A short examination of 
the formulas involved reveals inconsistencies with this approach, and again it is not 
surprising. We have been dealing with nonrelativistic quantum mechanical amplitudes, 
and they of course do not transform as scalars, but rather their squares transform as 
densities. Hence one must recast the momentum component expansion to exhibit this 
property: 

1 r rlk 
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In our example of Gaussian wave packets in the lab frame, this leads to 

$(^ S ) = (^fcW 2 - m2l »l (4) 

where the 4-vector a = (x Q , 0) contains the width x a of the packet in its rest frame. In 
a moving frame, one transforms to a' = ( , yx ,'yf3x ) which for nonrelativistic momenta 
does in fact lead to the expected Lorentz contraction of the wave packet. However, the 
general case leads to a much richer structure. 



The probability amplitude in the bound state rest frame becomes 
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where we have used initial lab wave packet momenta p a and pb, wave packet widths x Q 
and bound state width xb in their respective rest frames, and 
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are the approprate transformation factors. The integrals can be performed in the case 
of a nonrelativistic bound state ijixb 1, 

. , l2 const ( ; 2<Ti> 2 ) 

k,(t) 2 = r e V ^+(-*)^ 7 

where 

a = (mx B ) 2 + 2(mx ) 2 (l + 2s 2 1 + 2s 2 2 ) () 
b = (cisx - c 2 s 2 )(mxo) 2 

One finds that the effective separation time t s ® M = t s (P), i.e. the bound state "re- 
members" the momenta of the quark pair which led to its formation. For the situation 
Pa = Pb (where the hard production amplitude is maximum), one finds 



t QM = ^l [{mXB f + ( mXo ) 2 (2 + 4)] (0) 



Again, this is the effective separation time in the rest frame of the bound state, 
and the P-dependence is in addition to that which will occur in the Lorentz dilation 
transformation to the lab frame. 

Note that the direction of this factor is to increase the separation time as a func- 
tion of bound state momentum, i.e. to bring the quantum-mechanical parameters back 
into a region which could be compatible with experiment. For example, the parameters 
in Eq. 9 would lead to a reduction in P c from 5 GeV/c down to less than 3 GeV/c. 
No fit to the data is attempted here, since we are dealing only with a one-dimensional 
model. 

As mentioned previously, one must average over production position in a realistic 
nuclear geometry to get the final results. If we simulate this situation here with a uni- 
form density over a transverse size which yields the above P c 's, one finds a considerable 
flattening of the suppression curves. This effect tends to oppose the desirable results 
of the separation time increase, but in principle the total result will still be reflected 
in the data. 

An opposite point of view may also be examined in this context. There is a 
possibility that initial state effects in nuclear matter can mimic the P t dependence of 



suppression by skewing the transverse momentum distributions of the incoming partons 
which participate in the hard collisions. This scenario has been examined^ and under 
certain assumptions could possibly account for the entire effect observed in the NA38 
data. If this is the actual situation, one can use the results we have developed to put 
constraints on the plasma parameters such that it will not induce P t dependence in 
excess of that in the data. In our simple model, this can easily be seen as a requirement 
on the plasma lifetime t p > y/2t s , and also on the transverse size of the plasma d^$> t s . 

5. Summary 

The primary result of this investigation is that production of heavy quarkonium 
states in hadronic interactions must involve momentum components for the initially- 
produced quarks which are both relativistic and have an uncertainty which precludes 
using classical trajectories for their description. When applied to scenarios in which 
color forces are temporarily screened, such as in a quark gluon-plasma, one is led 
to a quantum mechanical version of a "formation time" (in reality it is a separation 
time), which must be calculated from overlap integrals of quantum mechanical wave 
packets and bound states. This has a significant effect on the Pt - dependence of 
suppression of quarkonium formation as interpreted as a signal for quark-gluon plasma 
formation. One unanticipated feature is the effect of the initial hard production process 
momentum dependence on the separation time, even when calculated in the bound state 
rest frame. It is shown that this additional effect tends to sharpen the cut-off slope of 
the suppression vs. P t curve. Averaging over nuclear geometry in the production process 
has a tendency to reduce its magnitude. Some remnants of this effect will be present 
in any scenario of plasma formation, even if initial state nuclear effects are responsible 
for the presently-observed Pt dependence in the data. Thus one can derive constraints 
on plasma parameters from either the observation or non-observation of this effect. 
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